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Построена самосогласованная со схемой возбуждения модель разрядного контура лазера на парах металла. Исследуются четы- 
ре электрические схемы накачки. Выполнен анализ влияния параметров разрядного контура на характеристики импульса накач- 
ки в лазере на парах свинца. 


1. Моделирование схем возбуждения 

На сегодняшний день объем продаж коммерчес- 
ких лазеров на парах металлов (ЛПМ) занимает нез- 
начительный сегмент общемирового рынка лазер- 
ных систем [1]. Под коммерческими, в данном слу- 
чае, подразумеваются лазеры с мощностью излуче- 
ния до 120 Вт и с КПД от 0,5 до 1 %. Такая ситуация 
сложилась в силу того, что основные производства 
ЛПМ ориентированы либо на научные исследова- 
ния, либо на создание мощных стационарных про- 
мышленных установок. Поэтому разработка ком- 
мерческих ЛПМ, является актуальной задачей. 

Поскольку ЛПМ еще не достаточно исследова- 
ны, а построение их действующих макетов связан- 
но со значительными финансовыми и временными 
затратами, то наиболее целесообразным является 
исследование данных лазеров методами численно- 
го моделирования. 
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Рис. 1. Исследуемые схемы разрядного контура 


Цель данной работы заключается в моделирова- 
нии характеристик разряда и параметров плазмы в 
течение импульса возбуждения в ЛПМ с использо- 
ванием самосогласованной нестационарной кине- 
тической модели. Следует отметить, что приведен- 
ная ниже модель разрядного контура справедлива 
для всех типов ЛПМ. Для того чтобы проводить мо- 
делирование необходимо задаться определенным 
рабочим веществом. В нашей работе таким вещест- 
вом выбраны пары свинца. Такой выбор обусловлен 
тем, что данному типу ЛПМ уделяется незначитель- 
ное внимание. По нашему мнению, исследование и 
оптимизация параметров лазера на парах свинца 
(ЛПС), на основе разработанной модели, позволить 
расширить области его применения. Математичес- 
кая модель ЛПС реализована в среде МаДаЪ 6.0. 


Исследованы основные схемы разрядного кон- 
тура ЛПС с: 

1. прямым разрядом накопительного конденсато- 
ра (рис. 1, а); 

2. шунтирующей индуктивностью (рис. 1 ,Ь)\ 

3. обостряющим конденсатором (рис. 1, с); 

4. шунтирующей индуктивностью и обостряю- 
щим конденсатором (рис. 1, сі). 


2. Математическая модель разрядного контура ЛПС 

Математическую модель ЛПС, состоящую из сис- 
темы жестких дифференциальных уравнений, приве- 
денных ниже, можно разделить на два связанных бло- 
ка, описывающих параметры электрической цепи 
возбуждения и параметры плазмы газоразрядной 
трубки (ГРТ). Уравнения первого блока, выражающие 
первый и второй законы Кирхгофа, будут отличаться 
для каждой из исследуемых схем. Уравнения второго 
блока для концентрации и температуры электронов 
будут одинаковыми для всех четырех схем. 

Уравнения Кирхгофа для разрядных контуров. 

1. Схема накачки с прямым разрядом накопи- 
тельного конденсатора Сі. 
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где Ц а - напряжение на накопительном конденса- 
торе Сі; /-ток разряда через ГРТ; к/ /Ѵ и Т е ) -актив- 
ное сопротивление разряда ГРТ, зависящее от кон- 
центрации Ы е и температуры Т е электронов; 
Л ы - индуктивность, представляющая собой сумму 
индуктивности разрядного контура Л к и индуктив- 
ности разряда С; Д„,(7) - динамическое сопротив- 
ление тиратрона, определяемое как: 


К, И (0 = КіНоІГ ех Р 
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где К, ко(/ и К іНо „ - сопротивления тиратрона в закры- 
том и открытом состояниях, С - время коммутации 
тиратрона. 

2. Схема накачки с шунтирующей индуктив- 
ностью. 
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где у=(Ь ! +Ь к )-( 1 1 1 +Ц)-Ц; І а - ток, протекающий че- 
рез накопительный конденсатор; /-ток, протекаю- 
щий через ГРТ; Ь„ - шунтирующая индуктивность. 
3. Схема накачки с обостряющим конденсатором. 
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- мощность, затрачиваемая на ионизацию атомов 
свинца; /, РЬ = 7,42 эВ - энергия ионизации атома 
свинца; 
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- мощность, затрачиваемая на нагрев тяжелых час- 
тиц в упругих соударениях; т е - масса электрона; 
т ѣ - масса атома свинца; м № - масса атома неона; 
- температура газа (равна температуре стенки ГРТ); 


9 -КГ 9 см V 


- скорость упругих столкновений электронов с 
атомами неона; 
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3 Т: 


- скорость упругих столкновений электронов с ио- 
нами; в - заряд электрона [3]; 
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где І Ср - ток, протекающий через обостряющий 
конденсатор. 

Уравнения, описывающие плазму ГРТ в момент 
импульса возбуждения. 

1. В данной модели использовано упрощенное 
выражение для концентрации электронов 7Ѵ„ кото- 
рое имеет вид: 
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- кулоновский логарифм; 
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- мощность, выделяющаяся при джоулевом нагре- 
ве;/ - плотность тока; р - удельное сопротивление 
плазмы; 
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где /Ѵрь - концентрация атомов свинца; 


к,-ръ = {° г и) ау т- р 


У Е* л 


-скорость ионизации атомов свинца, которая приня- 
та равной скорости возбуждения резонансных состоя- 
ний, что справедливо в режиме квазистационарной 
ионизации, когда каждый акт возбуждения сопровож- 
дается актом ионизации возбужденного состояния [2]; 
ѵ - скорость электронов; а г - сечение возбуждения ре- 
зонансного уровня свинца; ѵ, РЬ - множитель, характе- 
ризующий скорость ионизации свинца; /р Ь =4,37 эВ - 
энергия возбуждения резонансного уровня; 



- подгоночная формула, которая описывает зави- 
симость скорости возбуждения резонансного уров- 
ня свинца отГДЗ]. 


2. Выражение для температуры электронов, оп- 
ределяемое на основе уравнения теплового балан- 
са, имеет вид [3]: 


где 


Лі У. сІі 3 у 

й = ^іРЪ ' ^;РЬ ‘ ^РЬ 


- сопротивление плазмы активной среды лазера 
[4], где / - длина активной зоны ГРТ, У - площадь 
поперечного сечения активного элемента. С уче- 
том вклада кулоновских электрон-ионных и элект- 
рон-атомных столкновений сопротивление плаз- 
мы активной среды в предположении максвелловс- 
кой функции распределения электронов по энер- 
гиям можно записать в следующем виде [5]: 



где а ь а 2 - множители, определяющие вклад в соп- 
ротивление плазмы от электрон-атомных и элект- 
рон-ионных столкновений; г - радиус ГРТ. 


3. Результаты моделирования 

Используя разработанную модель разрядного 
контура ЛПС, проведен ряд исследований. Прежде 
всего, это сравнение полученных характеристик с 
экспериментальными, что позволяет сделать вывод 
об адекватности, построенной модели. Второе ис- 
следование заключалось в оценке влияния пара- 
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метров разрядного контура на характеристику им- 
пульса накачки. Данное исследование интересно 
потому, что именно параметры импульса накачки 
определяют степень инверсии населенностей ла- 
зерных уровней в ЛПС, а, следовательно, и мощ- 
ность выходного лазерного излучения. 

3.1. Сравнительный анализ результатов моделирования с 
экспериментом 

На рис. 2 -4 представлены характеристики раз- 
ряда и плазмы в течение импульса возбуждения в 
ЛПС, полученные в результате моделирования, ко- 
торые сравнивались с экспериментальными зави- 
симостями, заимствованными из работы [1]. Срав- 
нение проведено для схемы с прямым разрядом на- 
копительного конденсатора (рис. 1, а) 


I, А У, кВ 



Рис. 2. Временные зависимости напряжения и тока через 
ГРТ, для схемы 1; а), Ь) токи, с), сі) напряжения, рас- 
считанные нами и заимствованные из работы [1] со- 
ответственно 


Ие я ІО 1 ’ см ' 3 



Рис. 3. Концентрация электронов для схемы 1: а) расчет, Ь) 
эксперимент [1] 


Полученное расхождение результатов данной 
работы и приведенных в [ 1 ] объясняется использо- 
ванием в расчетах последних данных по констан- 
там скоростей процессов, а также построением бо- 
лее адекватной кинетической модели. 


Те. эВ 



Рис. 4. Изменение температуры электронов в процессе раз- 
ряда для схемы 1: а) расчет, Ь) эксперимент [1] 

3.2. Влияние параметров схемы возбуждения 

на импульс накачки ЛПС 

Как отмечалось выше, мощность выходного ла- 
зерного излучения определяется параметрами им- 
пульса накачки. 

3.2 Л. Влияние индуктивности. Показано, что 
при увеличении индуктивности контура возбужде- 
ния лазера на парах свинца амплитуда импульса то- 
ка через ГРТ уменьшается, одновременно с этим 
происходит увеличение его длительности. Таким 
образом, скорость нарастания импульса тока 
уменьшается. Соответственно уменьшается ско- 
рость нарастания напряженности электрического 
поля в разряде, а, следовательно, уменьшается и 
скорость нарастания температуры электронов. 

«.см 3 С ' 1 



Рис. 5. Зависимость константы скорости возбуждения от 
температуры электронов для: а) резонансного и Ь) 
метастабильного уровней атома свинца (линия излу- 
чения 722,9 нм) 

Сравнивая константы скоростей возбуждения 
резонансного (6р7в Т\) и метастабильного (6р 2 *Д) 
уровней, приведенных на рис. 3, видно, что при 
Т е <\,\ эВ эффективнее возбуждается метастабиль- 
ный уровень (6р 2 1 Д) атома свинца, а при превыше- 
нии этого значения эффективнее идет процесс воз- 
буждения резонансного уровня (6р7к Ч\). Назовем 
величину температуры электронов равную 1,1 эВ — 
пороговой. Следует отметить, что для эффективно- 
го возбуждения резонансного уровня необходимо, 
чтобы время, за которое Т е достигает порогового 
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значения, было много меньше некоторого крити- 
ческого значения і т определяемого из выражения: 

I I I 

\М 0 (1)М е (1)а г (0Ж=\М 0 (і№ е (!)а т (1)<11+\Л т ^ 

0 0 0 


где Ж 0 - концентрация атомов свинца в основном 
состоянии; а„ а т - скорости возбуждения резонанс- 
ного (г) и метастабильного ( т ) уровней; А т - вероят- 
ность спонтанного перехода с г-уровня на /я-уровень. 

Данное соотношение приближенно описывает 
процессы населения верхнего и нижнего лазерных 
уровней. Если за время, меньшее і а , температура 
электронов достигнет порогового значения, то в 
активной среде возникнет инверсия населенностей 
с последующей генерацией лазерного излучения, в 
противном случае генерация будет невозможна. 

Поскольку с увеличением индуктивности умень- 
шается скорость нарастания температуры электронов, 
то и время, за которое данная характеристика достиг- 
нет порогового значения, будет увеличиваться, что 
приведет к уменьшению степени инверсии населен- 
ностей уровней. В итоге средняя мощность лазерного 
излучения с увеличением индуктивности будет падать. 

Как было отмечено выше, кроме увеличения 
длительности переднего фронта импульса тока, уве- 
личение индуктивности приводит к увеличению его 
длительности. Увеличение времени импульса тока 
обеспечивает такой режим разряда, при котором воз- 
растает мощность, вкладываемая в плазму после ге- 
нерации лазерного излучения. Это приводит к уве- 
личению предымпульсных концентраций как элект- 
ронов, так и атомов свинца в метастабильном состо- 
янии. Соответственно, для достижения прежних 
предымпульсных значений концентраций электро- 
нов и атомов в метастабильном состоянии, необхо- 
димо увеличить время релаксации. Таким образом, 
увеличение индуктивности приводит к увеличению 
длительности импульса возбуждения и, соответ- 
ственно, времени релаксации, что означает умень- 
шение частоты следования импульсов. В свою оче- 
редь, ограничение частоты следования импульсов 
при неизменной энергии, вкладываемой в разряд за 
время импульса, приводит к ограничению средней 
энергии, вкладываемой в разряд, и соответственно к 
уменьшению мощности лазерного излучения. 

Следовательно, увеличение индуктивности, в 
конечном счете, приводит к уменьшению мощнос- 
ти лазерного излучения. Поэтому необходимо 
стремиться обеспечить как можно более низкое 
значение индуктивности разрядного контура. 

С другой стороны, величина индуктивности не 
должна быть меньше некоторого минимального зна- 
чения. Одним из требований, предъявляемых к пара- 
метрам импульса возбуждения, является условие апе- 
риодической разрядки накопительного конденсатора 
во время формирования импульса возбуждения, ко- 
торое выражается следующим неравенством [6]: 
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где Ь ІА - эквивалентная индуктивность разрядного 
контура, которая для схемы 1 является корректиру- 
ющей С к ; К іл - эквивалентное активное сопротив- 
ление разрядного контура; 

■^ЗсЬ =К(0 + *г', 

где Я л (1) - сопротивления тиратрона, і? раг - паразит- 
ное сопротивление проводов. 

Вышеприведенное неравенство определяет ус- 
ловие, накладываемое на минимальное значение 
величины корректирующей индуктивности: 

4 >4(Ль(0 + ^ Р аг) 2 С 8 - 

3.2.2. Влияние емкости накопительного конденса- 
тора. Показано, что при увеличении емкости нако- 
пительного конденсатора увеличивается длитель- 
ность и амплитуда импульса тока. Крутизна перед- 
него фронта импульса практически не изменяется. 

Увеличение длительности импульса тока при- 
водит к необходимости увеличения времени релак- 
сации, а, следовательно, и ограничению частоты 
следования импульсов возбуждения, что приводит 
к уменьшению мощности лазерного излучения. 
Механизм ограничения частоты импульсов воз- 
буждения и уменьшения лазерного излучения рас- 
смотрен в предыдущем пункте. В настоящее время 
ведутся работы по оптимизации лазеров на парах 
металлов, в частности лазеров на парах меди и 
свинца, для достижения максимально возможной 
частоты повторения импульсов. Увеличение часто- 
ты следования импульсов, сопровождающееся 
уменьшением мощности лазерного излучения, при 
увеличении емкости накопительного конденсато- 
ра, определяет максимальное значение последнего. 

С другой стороны величина емкости определяет 
энергию, которая запасается в накопительном кон- 
денсаторе и часть которой вкладывается в разряд. 
Уменьшение величины емкости приводит к уменьше- 
нию энергии, запасаемой в накопительном конденса- 
торе, а, следовательно, и энергии, вкладываемой в 
разряд. Данное условие определяет минимальную ве- 
личину емкости накопительного конденсатора. 

3.2.3. Влияние величины начального напряжения 
на накопительном конденсаторе. Показано, что уве- 
личение начального напряжения на накопитель- 
ном конденсаторе увеличивает амплитуду и кру- 
тизну переднего фронта импульса тока через ГРТ 
при неизменной длительности. С учетом предыду- 
щих пунктов следует, что начальное напряжение на 
накопительном конденсаторе должно быть как 
можно большим. Ограничение на величину на- 
чального напряжения обусловлено свойствами 
коммутирующего элемента: максимально допусти- 
мой коммутируемой мощности и скорости нарас- 
тания импульса тока и напряжения. Превышение 
этих значений недопустимо, так как приводит к 
уменьшению срока службы и выходу из строя ком- 
мутирующего элемента. Поэтому необходимо вы- 
бирать коммутирующие устройства с высокими па- 
раметрами максимально допустимой коммутируе- 
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мой мощности, скорости нарастания тока и напря- 
жения. Наилучшие параметры имеют мощные по- 
лупроводниковые коммутаторы, такие как тирис- 
торы, ІСВТ и М08РЕТ транзисторы. 

Сравнивая блоки возбуждения, построенные с 
использованием вакуумных и полупроводниковых 
коммутаторов, следует отметить, что последние об- 
ладают более высоким сроком службы, большим 
КПД, а также лучшими условиями синхронизации 
коммутирующих элементов при построении слож- 
ных систем [7]. 

3.2.4. Влияние начальной концентрации электро- 
нов. Показано, что при увеличении начальной кон- 
центрации электронов, сопротивление плазмы 
уменьшается, т.е. повышается проводимость кана- 
ла. Это препятствует эффективному вводу энергии 
в разряд и, соответственно, снижает максимально 
достижимую в течение импульса накачки темпера- 
туру электронов. 

Механизм влияния предымпульсной концент- 
рации электронов на процессы, протекающие в те- 
чение импульса возбуждения можно, описать сле- 
дующим образом. Концентрация электронов опре- 
деляет теплоемкость электронного газа, температу- 
ру которого необходимо повысить до определенно- 
го значения, определяемого преобладанием воз- 
буждения резонансных уровней по сравнению с 
метастабильными (см. выше). При наличии реаль- 
ного источника питания, выходная мощность ко- 
торого является ограниченной величиной, превы- 
шение определенной предымпульсной концентра- 
ции электронов приведет к невозможности реали- 
зации достаточного разогрева электронного газа и, 
как следствие, возможности существования инвер- 
сии в активной среде лазера. Таким образом, необ- 
ходимое условие возникновения лазерного излуче- 
ния заключается в том, что начальная концентра- 
ция электронов должна быть меньше некоторого 
критического значения [3]. 

3.2.5. Влияние шунтирующей индуктивности. Для 
схемы 2, ток через ГРТ практически не зависит от 
величины шунтирующей индуктивности. 

Определим величину шунтирующей индуктив- 
ности. Исходя из приведенного в пп. 3.2.1 условия 
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апериодического разряда накопительного конден- 
сатора, а также с учетом того, что для схемы 2 
Ад » 4 + А [6], получаем: 

А ~ 4 ( ( 0 + я раг ) 2 С з ~ Ь к > 

где А - шунтирующая индуктивность. 

3.2.6. Влияние емкости обостряющего конденса- 
тора. Для схемы 3, при увеличении емкости обост- 
ряющего конденсатора скорость нарастания перед- 
него фронта тока падает. Это объясняется тем, что 
с увеличением емкости увеличивается постоянная 
времени цепи разряда. 

Уменьшение крутизны переднего фронта им- 
пульса, как показано выше (см. п. 3.2.1), является 
нежелательным, поскольку это приводит к умень- 
шению мощности выходного лазерного излучения. 

3.3. Влияние типа схемы возбуждения 

Сравнивая полученные характеристики разряда 
и плазмы четырех схем между собой можно прийти 
к выводу, что наиболее предпочтительными схема- 
ми накачки, являются схемы 3 и 4, т.к. в этом слу- 
чае увеличиваются скорости нарастания и макси- 
мальные значения импульса тока по сравнению с 
схемами 1 и 2. 

Заключение 

1 . Построена самосогласованная со схемой возбуж- 
дения модель лазера на парах металла. Модель 
включает в себя уравнения Кирхгофа, описываю- 
щие схемы накачки, и уравнения для определе- 
ния температуры и концентрации электронов. 

2. Исследованы четыре электрические схемы на- 
качки. Выполнен анализ влияния параметров 
разрядных контуров на характеристики им- 
пульса накачки. 

3. Разработанные математическая модель ЛПМ и 
программа расчета основных параметров рабо- 
ты лазера на парах свинца позволяют произво- 
дить исследование ЛПМ, а также оптимизацию 
параметров устройства в целом и сократить фи- 
нансовые и временные затраты при разработке 
новых технологических образцов ЛПМ. 
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